4D Spin Glasses in Magnetic Field Have a Mean Field like Phase 
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By using numerical simulations we show that the AD J = ±1 Edwards Anderson spin glass in 
magnetic field undergoes a mean field like phase transition. We use a dynamical approach: we 
simulate large lattices (of volume V) and work out the behavior of the system in limit where both 
t and V go to infinity, but where the limit V — > oo is taken first. By showing that the dynamic 
overlap q converges to a value smaller than the static one we exhibit replica symmetry breaking. 
The critical exponents are compatible with the ones obtained by mean field computations. 



The mean field solution of spin glass systems Q| con- 
tains many new features. It tells us that systems with 
quenched disorder can have a large number of stable 
states, not related by explicit symmetries of the origi- 
nal Hamiltonian, and that the space of these states is 
embedded with an ultrametric structure. Moreover, the 
system stays critical for all T < T c and the phase tran- 
sition of Replica Symmetry Breaking (RSB) survives the 
presence of a finite magnetic field h. 

The mean field paradigm needs to be analyzed, in or- 
der to understand how many of its peculiar features are 
shared by the finite dimensional, physically relevant case. 
In spite of the technical difficulties, in the last years many 
progresses have been done. It is for example remarkable 
that recent rigorous results seem to support strongly 
(after some initial different feelings ||) the viability of 
the mean field approach for the description of finite di- 
mensional systems. It has been shown H that the rig- 
orous finite dimensional construction of |3| leads to self- 
averaging quantities exactly where the mean field con- 
struction would also produce self-averaging observables, 
and Guerra has shown that the main part (and maybe 
all) of the replica predictions on the fluctuations of non- 
self-averaging quantities applies to the broken phase of 
finite dimensional disordered systems (see also Q). 

Monte Carlo simulations are an important tool to es- 
tablish how much of the mean field description survives 
in the finite dimensional case |5|||. For example there 
is now evidence for the existence of a 3D mean-field like 
critical point, for the existence of an ultrametric struc- 
ture in 4D, and for a dynamical behavior of finite di- 
mensional systems very similar to the one that can be 
found analytically in the Sherrington-Kirkpatrick mean 
field model. 

The question of the existence of a de Almeida- Thou- 
less M line, i.e. of the existence of a phase transition 



in finite magnetic field, is maybe the most relevant open 
problem. Even if a large amount of numerical work has 
been done to clarify this issue , a clear cut answer is 
still lacking. Most of the numerical work suggests that a 
transition exists (even if some studies suggest the oppo- 
site conclusion), but the question is a very delicate one: 
one finds probability distributions that do not have a 
very clear behavior, and it is very difficult to thermalize 
large systems in the low temperature (T) region. Even 
the most recent numerical work of || does not reach un- 
ambiguous conclusions. 

Here we hope to settle the question, by showing in a 
non-ambiguous way that the 4D spin glass with quenched 
couplings J = ±1 in finite magnetic field undergoes a 
mean field like phase transition. 

We use a dynamical approach. If a large system is 
cooled down to a temperature Tf, starting from the high 
temperature region, after a time t the correlation func- 
tions are different from zero (in a statistically significant 
way) only up to distances smaller than a dynamic corre- 
lation length £(t). Often (and this seems to be the case 
of spin glasses in the low T phase) increases as a 
power of t, i.e. £(t) oc t 1 ^^. If the lattice size L is 
larger than for large times the system is locally, but 
not globally thermalized: in the case of an infinite lattice 
this is always the case, independently from the value of 
t. As we shall see later our choice of the lattice volume, 
V — 20 4 , is such that we stay in this situation. Then 
by using power fits we keep the large time limit under 
control, and we determine with high precision the infi- 
nite time expectation values, always in the phase where 

A key prediction of the theory of replica symmetry 
breaking is that when comparing different realizations of 
the system we find that there are local quantities which 
take a different value in the t — > oo and V — ► oo limits in 
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the two regions L ^> and L <C £(t). In the following 
we will call respectively dynamic and static the expecta- 
tion values computed in the first and in the second region. 
We aim to show that in four dimensional spin glasses in 
magnetic field at low temperature the two expectation 
values are different and therefore the replica symmetry is 
broken, as expected from mean field computations. 

This work contains two kind of results. First of all 
we discuss some inequalities, both at finite ft and in the 
h — > limit, that can be violated only if replica symmetry 
is broken. We use our numerical simulations to show that 
indeed such inequalities are broken for T small enough. 
Second we show that our data for the overlap and for the 
underlying time scales obey an impressive scaling versus 
the magnetic field, and that the critical exponents turn 
out to be very similar to the mean field theoretical pre- 
diction. 

Numerical data are drawn from dynamical runs sched- 
uled according to the following scheme. We start at 
T\ = 3.00 (the value of the critical temperature at h = 
is close to 2.0) and decrease T with steps of 0.25 down to 
Tg = 1.00. For an annealing run of level k at T — T n , n-2 k 
steps are performed with k spanning from 5 to 16 (14 for 
the larger fields) and ft from 0.2 to 0.6 (runs at h = 0.1 do 
not reach enough precision to be used for fitting, and have 
only been included in the matching analysis, see later). 
In total our longer runs involve order of 3 millions sweeps 
of the 20 4 lattice. The quantity t = 2 k plays the role of a 
time: we will be extrapolating expectation values on k at 
given T and ft. We use two copies of the system a and r 
in each realization of the quenched couplings to compute 
the overlap q = ^jCjTj/V. We use a multispin coded 
algorithm jl0| that allows to flip more than 70 ■ 10 6 spins 
per second on a Digital a workstation 500/333. We have 
averaged over 64 samples for each k and ft value (for a 
very few cases we only have 32 samples). In principle we 
could also put the system at the final temperature by a 
a sudden quench. We have followed the previous proce- 
dure for two reasons: (a) Finite time effects are smaller 
and the infinite time extrapolation is easier, (b) We can 
collect in one run data at different temperatures. 

The first kind of evidence is based on our results for 
q(t) at fixed ft. We extrapolate q(t) to its value for infinite 
time q D (q D , where D stands for dynamical, is also the 
minimum allowed value for q at equilibrium, q mm |H). 
We compare q D with the static value q computed with 
equilibrium runs for a 7 4 system B (preliminary runs 
of |ll]] confirm the determination of q s of |J down to 

r = i.oo). 

We show that at low T q D < q s strictly, i.e. that 
replica symmetry is broken. In fig. (Q) we show two 
typical fits for low T = 1.0, at h = 0.4 and h = 0.2 (here 
we are using for fitting all the data points: see later for 
scaling time windows). Moreover, since the values of q 
increase with the lattice sizes in the static runs of Hj , our 
evidence is safe also from the point of view of finite size 



effects. 




10 J 

Annealing time 
FIG. 1. q(t) versus t and best fit. 

The power fits are very good. One finds qif) = 0.56(1) 
+ 0.44(3) t- 28 ( 3 ) and q(t) = 0.50(2) + 0.58(2) t-- 21 ( 3 ) 
respectively at ft, = 0.4 and 0.2. Both fits have a very 
good x 2 - 

In figure (j|) we plot our data for q( D ) (dashed curve 
and error bars) at ft, = 0.4 and the static data of ||. For 
high values of T the data are in perfect agreement, while 
at T = 1.5 the two curves start to split in a statistically 
significant way. 
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2. q D from our runs 


and q s from M versus T. 



Both at T = 1.25 and at T = 1.0 it is clear that the 
difference of the dynamic and the static value is both 
statistically and systematically significant with a large 
confidence level. So at ft = 0.4 we have evidence that for 
T < 1.25 the system is in a mean field like broken phase. 
These results are in good agreement with the data of || 
that suggests a transition near T — 1.5 at this value of 
the magnetic field. 
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The reader could wonder if in our simulations the in- 
equality L » £(t) is satisfied. By estimating the expo- 
nent z(T) from the simulations at h = we find that 
z » 10 at T = 1.0 and that the bound should be satu- 
rated for times O(10 13 ), which is much larger than the 
largest times scales of our numerical simulation. Simple 
power fits to energy, overlap and magnetization and to 
their fluctuations are good: this is an independent indi- 
cation of the fact that that all our data are in the region 
where L>^(t). Moreover, even if L was close to our 
conclusion would be strengthened since then the differ- 
ence among the measured dynamical value and the static 
value could only decrease. 

Next we discuss the h — > limit. We consider the sus- 
ceptibility x = linih^o X ■ ^ replica symmetry is realized 
we have that for h — > the susceptibility \ = (3(1 — q), 
where there is no ambiguity in the definition of q. We 
can thus define 



P(0)/2 (the value of P(0) at this temperature is about 
•5 M). 



q(h) 



1 — T 



m(h) 



(1) 



If in the h — > limit q ^ q D , than replica symmetry 
is broken. In a theory where replica symmetry is broken 
the small h limit of q is q s . More precisely for finite h 
we find that q(h) = q s (h) + 0(h 2 ). In fig.(|].a) we show 
q D at T = 1.0 as a function of h (empty dots), together 
with the values of q{h) (filled dots). The two functions 
do not extrapolate to the same value at h = 0. Replica 
symmetry is broken in the region where q D is smaller 
than the h = limit of q. If we neglect terms of order h 2 
(the difference among q and q s is about .02 at h = .4), 
replica symmetry must be broken when the two curves 
differs in a statistically significant way (i.e., in our case 
at T = 1.0, for h < 0.5). 

As we will discuss later we have determined a rescaling 
of times as a function of h that makes the curves q(t) at 
different fields universal. We can thus determine consis- 
tent h dependent time windows, that allow us to com- 
pare homogeneous time regimes at the different h values. 
The fc-windows used for these scaling fits (t = 2 k ) are 
8 - 16 at h = 0.2, 7 - 14 at h = 0.3, 6 - 13 at h = 0.4, 
5 - 11 at h = 0.5 and 5 - 10 at h = 0.6. The results of 
the corresponding fits are given in figure (@.b). Here the 
points have somehow a larger error (since we use less data 
point for fitting) but we expect the systematic error to be 
smaller. The points at h — 0.2, for example, appear more 
consistent thanks to the elimination of short time effects. 
The emerging physical picture is independent from the 
fitting scheme. 

For small h in the SK model q(h) — q(0) + Rq 
R ~ P(0)/2, and P(0) is the value of the function P(q) 
at q = when h = 0. We have found that in the region 
where q ^ q D the data are compatible with a quadratical 
dependence over q D . We find for example R = .4 at 
T = 1.25, which is of the same order of magnitude of 



I,, where 
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(b): T=1.0, rescaled windows 
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FIG. 3. q D and q from fits including all time points and 
from fits done over rescaled time windows. 

Let us give some information about the exponent of the 
power-law fit we have determined for the decay of q{t). 
At low T (1.0 and 1.25) such exponents are between .2 
and .3 with no apparent systematic dependence on the 
magnetic field. The fits on rescaled time windows give 
higher values than the fits on all points (basically for low 
T values the results are fixed around We have also 
fitted the energy with a power decay to its asymptotic 
value. Here the decay exponent can be estimated with 
good precision, and for low T it does not depend on h. 
For example at T = 1.00 we find an exponent of .435 for 
h going from 0.1 to 0.4. At T — 1.25 it is 0.47 for h going 
from 0.1 to 0.3, while at T — 1.5 it already has a small 
dependence on h. 

Our last numerical evidence is based on rescaling the 
functions quit) obtained for different values of h. We have 
rescaled our data obtained at different magnetic field val- 
ues according to qh'it') = A(h)qh(B(h)t). 

The coefficients for the rescaling to a fixed value of h! 
are fitted to the form A(h) ~ Ah", and B(h) ~ Bh T . 

The value of the crossover exponents may be found 
by assuming that the quantity J d x h(x) 2 q(x) is dimen- 
sionless: the coupling term in replica space has a form 
/ d 4 x h 2 (x) J2 a b Qa,b(x)- The dimension of q in the dy- 
namic approach can be reconstructed by the decay of the 
correlation functions. An approximate formula (which 
seems to work reasonably well |l2|]) has been proposed 
in H|: (q x q ) cx x~ A , with A w (D - 2)/2. This formula, 
together with dimensional analysis, implies that in D = 4 
u> ~ j. A similar analysis shows that in D = 3 A = jj. 
Since the dynamical critical exponent is of order 8 in the 
temperature range around T = 1.25 the same argument 
implies that r ~ 4. 

We have determined the best coefficients A and B by 
minimizing the square difference of the two functions. 
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Below we discuss results at T = 1.0. The fits are very 
good: all the rescaled q have ratios systematically com- 
patible with 1, and there is no need for a further extrap- 
olation or corrections to our scaling formula. 
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FIG. 4. Rescaled functions q(t) versus annealing time 



In figure (Q) we show the rescaled functions (horizon- 
tal and vertical scales are given by the fact that we have 
kept fixed the values at h = 0.2). The scaling is obeyed 
remarkably well: it works over six time decades, and in a 
range of magnetic fields going from .1 to .6. The errors on 
data point are not plotted since they would blur the fig- 
ure. They are of the order of 0.01-0.02. For example the 
point at h = 0.1 with largest t values, that is slightly out 
of the enveloping curve, is statistically compatible with 
the other points. A and B determined with the fits of fig- 
ure can now be fitted with power laws. We plot them 
in figure (|^), with A represented by the upper points and 
B from the lower ones, together with the best fits (the 
fitting function are normalized in such a way to give 1 
at h = 0.6). The best fit gives A(h) ~ 1.14(l)/r 25 W, 
and B{h) ~ 5.4(3)/i 33(1) . Even if this is a qualitative 
test, since we have only a rough estimate from the mean 
field approach, and in this case we have not analyzed the 
statistical and systematic error in great detail (since sys- 
tematic error could be quite large for this measurement) 
the agreement with the values one would expect from the 
mean field solution turns out to be remarkably good. 

We acknowledge useful discussions with C. Naitza, F. 
Ricci-Tersenghi and J. J. Ruiz-Lorenzo. We warmly 
thank F. Ritort for interesting correspondence and com- 
ments. 
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